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The Jaynes-Cummings model, describing the interaction between a single two-level system and 
a photonic mode, has been used to describe a large variety of systems, ranging from cavity quantum 
electrodynamics, trapped ions, to superconducting qubits coupled to resonators. Recently there has been 
renewed interest in studying the quantum strong-coupling (QSC) regime, where states with photon number 
greater than one are excited. This regime has been recently achieved in semiconductor nanostructures, 
where a quantum dot is trapped in a planar microcavity. Here we study the quantum strong- coupling regime 
by calculating its photoluminescence (PL) properties under a pulsed excitation. We discuss the changes in 
the PL as the QSC regime is reached, which transitions between a peak around the cavity resonance to a 
doublet. We particularly examine the variations of the PL in the time domain, under regimes of short and 
long pulse times relative to the microcavity decay time. 

Strong coupling is a phenomenon that has gained fundamental importance in the field of quantum tech- 
nologies, where combinations of different systems are coupled quantum mechanically. Typically of such 
engineered systems, a matter field consisting of two levels is coupled to a photonic or phononic bosonic 
mode. For low densities, defined as being less than one excitation per two-level system, the new excitations of the 
system are superpositions of the matter and bosonic fields. The ability to create a superposition between matter 
and bosonic fields is important in terms of quantum information processing applications due to their ability to 
couple isolated quantum systems together. For example, a common bosonic mode may be used as a quantum 
communication channel between qubits, providing effective two-qubit interactions. The Jaynes-Cummings 
Hamiltonian, which describes any two-level system interacting with a single bosonic mode, is given by 



with 



Hmt=g 



[g a^ + G^a) 



(1) 

(2) 
(3) 



where g- are the raising and lowering operators for the two-level system, and is the creation operator for the 
bosonic mode. The coupling between the photon and the two level system is given by ^, the transition frequency of 
the two-level system is co, and A is the detuning between the bosonic mode and the qubit^"^. 

At high densities, where there is more than one excitation per two-level system, the model predicts the 
quantum strong- coupling (QSC) regime. The quantum strong- coupling regime has been studied using super- 
conducting qubits^"^. This may be seen by directly diagonalizing the Hamiltonian, which may be easily performed 
since it preserves the total excitation number tixox = ci'a + |X)(X| . Considering the case of a resonant (A = 0) qubit 
and cavity system, the eigenstates may be found to be 

\±,n)^±{\X}\n-l}±\G}\n}) (4) 

where \n) = —=(^a^y\0). The ground and excited states of the two-level system are denoted |G) and 
respectively. The energies of the states are 
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E^=nhw±2g^ (5) 

thus for each excitation number n (with the exception o{n = 0 which 
is unique) there are two solutions forming the Jaynes-Cummings 
ladder, with the energy splitting increasing with the excitation num- 
ber n. Recently an experiment studying the Jaynes-Cummings ladder 
in the quantum strong- coupling regime was realized, where a single 
quantum dot was placed in a microcavity structure^'^. The highly 
excited states of the Jaynes-Cummings ladder were accessed by a 
combination of off-resonant pumping which first excites the ^tot = 
1 manifold, then resonant pumping which results in climbing to 
higher /itot states. 

In this paper we perform a numerical study of the photolumines- 
cence of the quantum dot-microcavity system. While similar studies 
of the photoluminescence have been made for the Jaynes-Cummings 
model in a number of works^"^\ the pumping scheme that is used 
in these works typically consider a continuous -wave pumping. 
However, in order to achieve the high densities required for the 
QSC regime, experiments typically use a pulsed excitation due to 
the power constraints of the typical lasers used. Such pulsed excita- 
tions give a time -dependent, non- equilibrium element to the prob- 
lem. The dynamics of the system is determined not only by the 
Hamiltonian (1), but due to the open-dissipative nature of the matter 
and bosonic fields, which have a finite decay time. In our calculations, 
we consider both regimes where the pulse duration is shorter and 
longer than the lifetime of the particles in the system. Due to the 
pulsed excitation, the system is in general always in a non-equilib- 
rium state. This allows for the possibility of directly studying the 
system during the decay process via observation of the PL character- 
istics. 

We also consider the effects of various pumping configurations, 
namely with respect to either the matter or bosonic fields^^'^^. In 
general, it is also possible to observe the photoluminescence that is 
associated with either the matter or bosonic fields. In the case of 
experiments in Ref 7, either the PL due to the photons leaking out 
of the cavity, or the PL due to recombination of the excitons can be 
distinguished. These can in principle give different spectral charac- 
teristics, as we find explicitly from our calculations. This gives four 
possible pumping/measurement schemes, the relevant one depends 
on the particular experimental configuration used. We calculate all 
four cases and find that there are particular spectral characteristics 
common to all configurations which give spectral evidence for the 
QSC regime being reached. 

Results 

To model the open-dissipative dynamics of the coupled cavity- 
quantum dot system, we use a master equation with the coherent 

T3 

c 
3 




Figure 1 | The particles, oscillating between the ground state | G) and the 
excited state | X), are pumped in the cavity by a intense pump laser, either 

Pc or Px. The excited particles in the cavity decay with either the exciton 
decay rate y or cavity photon decay rate k. The coupling strength between 
the cavity photon and a two-level atom is given by g. 



part obeying (1), Lindblad photon and exciton decay terms with rates 
K and y respectively, and Lindblad photon and exciton pump terms 
with rates Pc and Px respectively (see Figure 1). The pump terms are 
assumed to obey a Gaussian profile in the time domain with time 
duration Tp (see Methods). The density matrix elements are numer- 
ically evolved in time, and the photoluminescence is calculated using 
two-time correlation functions evaluated using the quantum regres- 
sion theorem^"^\ A cutoff n^ax for the number of photons is made 
sufficiently large such that the trace of the density matrix is close to 
unity for all times. For concreteness we have chosen parameters that 
are of the order of the experiment performed in Ref 7. Specifically, 
the coupling constant^ is of the order of 10 meV, while the excitation 
energy fiw is of the order of eV; thus we fix g= 10 and fiw= 1000 
where all the parameters are in units of meV. The detuning para- 
meter is set to A = 0 throughout this study. 

Density matrix dynamics. We consider two regimes in this study, 
where the pump timescale Tp is either much longer or shorter than 
the polariton lifetime 1/k, when the system is pumped by either Pc or 
Px' Considering the case of long pumping (Tp > 1/k) first, we show 
the time evolution of the populations and the pump -pulse profile in 
Figure 2. The pumping pulse is centered at to = 0, and the system 
starts in the vacuum state \G)\n = 0). We see that as the pump is 
increased, the probability that the state is in the vacuum state 
decreases, resulting in excitations being created. Cutoffs of n^^^ = 
50 for Pc and /imax = 30 for Px are sufficient for the density matrix to 
have Trp = 1 for all times. At the strongest pumping intensity, there 
is a dip in some of the probabilities (e.g. the \X)\n = 1) state), 
signifying that excitations beyond the ^tot — 1 manifold are being 
excited. This is by definition the QSC regime, and confirms that with 
suitably intense pumping such regimes may be reached. 

Similarly, Figure 3 shows the density matrix elements for the case 
of a short pulse time Tp < 1/k. Again we see that many of the states 
above the ^tot — 1 manifold are being excited, confirming that either 
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Figure 2 | Density matrix dynamics /?{t) for long pulse pumping with Tp 
> \Ik. The time evolution of: (a) the occupation PGn,Gn for the ground 
state On under photon pumping Pq (b) the occupation Pxn,xn for the 
excited state Xn under photon pumping Pq (c) the occupation PGn,Gn for 
the ground state On under exciton pumping Px, and (d) the occupation 
Pxn,xn for the excited state Xn under exciton pumping P^. Curves are 
plotted in order of red {n = 0), green {n = 1), blue {n = 2), red (n = 3), 
green (/7 = 4), blue {n= 5), etc., with photon number n. As references, the 
dashed curves indicate the pumping profile. Parameters used here are g = 
10, ^ft)= 1000, Tp = 10, y = 0, and k = I. Negative times refer to prior to 
the pulse at t = 0. 
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Figure 3 | Density matrix dynamics /?{t) for short pumping pulse with Tp 
< \Ik. The time evolution of: (a) the occupation pGn,Gn for the ground 
state Gn under photon pumping Pq (b) the occupation Pxn,xn for the 
excited state Xn under photon pumping Pq (c) the occupation pGn,Gn for 
the ground state Gn under exciton pumping Px, and (d) the occupation 
Pxn,xn for the excited state Xn under exciton pumping P^. Curves are 
plotted in order of red {n = 0), green {n = 1), blue {n = 2), red {n = 3), 
green {n = 4), blue {n= 5), etc., with photon number n. As references, the 
dashed curves indicate the pumping profile. Unlike the case with long 
pulse, the system evolves asymmetrically with time. Parameters used here 
are^= 10, ?2Ct)= 1000, Tp= 1.5, y = 0, and k: = 0.1. Negative times refer to 
prior to the pulse at t = 0. 



Cavity photon 



pumping scheme may be used to reach the QSC regime. We however 
see an asymmetric behavior to the density matrix elements, unlike 
Figure 2 which was approximately symmetric. This can be explained 
due to the relative timescales of the pumping, where either the system 
can respond fast enough to changes in the pumping (the former 
case), and the pump is faster than the timescale of the system (the 
latter case). In the case of the short pulse, the polaritons require a 
time ~ 1/k; before they can escape the system, thus the decay process 
itself is being observed for times t > to. 

Photoluminescence spectra. We now show the photoluminescence 
dynamics of the quantum dot system. We consider two cases of long 
pump pulses (ip > 1/k, 1/y) and short pump pulses (ip < 1/k, 1/y). 
For each case we evaluate the four cases of cavity photon pumping 
and exciton pumping with the measurement of the photonic or 
excitonic PL (Figure 4 and 6). From the evolved matrix elements 
we also calculate the corresponding expectation values for the 
photon number = a^a and the exciton number = \X){X\ 
(Figure 5 and 7). We assume that the exciton decay rate y is 
negligible when measuring the photon mode, whereas the cavity 
photon decay rate k is negligible when measuring the exciton 
mode, such that y = 0 and k: = 0 for photonic PL and excitonic 
PL, respectively. 

Photon PL and photon pumping. First we examine the case where the 
system is photon pumped, and decays via photon emission k = 1. 
The average photon and exciton number is shown in Figure 5(a) for 
the long pump pulse case. We see that the photon number is much 
greater than the exciton number at the peak densities, which is a 
consequence of there being a maximum exciton number of 1, 
whereas there is no limit to the maximum photon number. For the 
quantum dot system this is a statement of Pauli's exclusion principle, 
as two excitons may not occupy the same dot (assuming the 
dot size is the same or smaller than the exciton Bohr radius). 
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Figure 4 | Photoluminescence spectra for long pump pulses. The time evolution of: (a) photonic PL by photon pumping Pq (b) photonic PL by exciton 
pumping Px, (c) excitonic PL by photon pumping Pq and (d) excitonic PL by exciton pumping Px in the case of long pump pulse. The peak PL 
approaches the cavity photon mode at high densities, whereas two peaks at UP and LP are brighter at low densities. Dashed line indicates the cavity photon 
energy. Parameters used are g = 10,hco= 1000, Tp= 10, and: (a) PJ^^ = 0.8, k:= l,y= 0,n^^^= 50; (h) P^^ = 6.0, k= l,y= 0, n^^^= 30; (c) P^^^ = 0.3, 
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Figure 5 | The mean number of cavity photons (n^h) and excitons (/igx) for long pump pulses. The cases shown are (a) photonic PL by photon pumping 
Pq (b) photonic PL by exciton pumping Px> (c) excitonic PL by photon pumping Pq and (d) excitonic PL by exciton pumping P^in the case of long 



pump pulse. Parameters used are ^ = 10, fico= 1000, Tp= 10, and: (a) P^^ 
(c) Pg^^ = 0.3, K= 0, y= 1, nmax= 120; and (d) PJ^^ = 3.5, k= 0, y= 1, n„ 

At high densities, there are therefore many more photons than exci- 
tons, again illustrating the QSC regime. For very high den- 
sities, such that (/iph) > 1, the bosonic operators of the 
Hamiltonian (1) may be written as a c-number, giving the effective 
Hamiltonian^^ 
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= 40. 
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Figure 6 | Photoluminescence spectra for short pump pulses. The time evolution of: (a) photonic PL by photon pumping Pq (b) photonic PL by exciton 
pumping Px, (c) excitonic PL by photon pumping Po and (d) excitonic PL by exciton pumping P^. The peak PL approaches to the cavity photon mode at 
high densities, whereas two peaks at UP and LP are bright in addition to the peaks near the cavity photon resonance at low densities. Dashed line indicates 
the cavity photon energy. Parameters used are ^ = 10, ?za>= 1000, Tp= 1.5, and: (a) Pc^"" = 0.45, k:= 0.1, y= 0, ^max= 40; (b) PJ'^'' = 1.3, k:= 0.1, y= 0, 
^max= 30; (c) P^^^ = 0.45, K= 0, y= 0.1, n^^^= 40; and (d) PJ^^ = 1.3, k= 0, y= 0.1, n^^^= 40. 
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Figure 7 | The mean number of cavity photons (n^h) and excitons (/igx) for short pump pulses. The cases shown are (a) photonic PL and photon 
pumping Po (b) photonic PL and exciton pumping Px> (c) excitonic PL and photon pumping Po and (d) excitonic PL and exciton pumping P^. 
Parameters used are^ = 10, hco= 1000, Tp= 1.5, and: (a) P^^ = 0A5, k= 0.1, y= 0, n^^^= 40; (b) PJ^^ = 1.3, k= 0.1, y= 0, n^^^= 30; (c) P^^^ = 0.45, 
K= 0, y= 0.1, nmax= 40; and (d) PJ^^ = 1.3, k:= 0, y= 0.1, ^n,ax= 40. 



see that the exciton number never exceeds this value, in agreement 
with this argument. 

In Figure 4(a) we see that at the highest densities the brightest peak 
comes from the cavity photon energy. As the system decays via 
leakage through the microcavity, the peak moves towards the lower 
polariton (LP) and upper polariton (UP) resonances, finally resulting 
in only the LP and UP resonances at the lowest density. A qualitative 
understanding of this spectrum can be described by using the dressed 
states \ ±, n) defined in equation (4). The strong photon pumping 
initially creates a state which is 

(at)"|G)|0)x|G)|n) 

1 (7) 

= —(\+,n)-\-,n)). 

As the photon PL is observed due to a loss of a photon out of the 
system, this state transitions to 

a\-\- ,n) — ^/n\-\- ,n — l), 

(8) 

a\-,n)^^/n\-,n-l), 

where we have assumed 1. The transition energy is then 

A£=£±-£±_, (9) 

which corresponds to the cavity photon energy. In the low density 
limit, we have energy peaks at UP and LP as shown in Figure 4(a). 
Thus, we find that the energy at the cavity photon arises from the 
enhanced transitions of| + ,n)^| + ,n — 1) and \ — ,n)^\ — ,n — l) 
due to the increased number of photons. The initial state has the 
same amplitude for the states | + , n) and | — , n), which is the reason 
for the symmetrical energy peaks to the cavity photon mode. 

In Figure 7 we plot the average number of photons and excitons for 
the short pump pulse case. We again see an asymmetric distribution 
of the particle numbers with respect to the maximum pump intens- 
ity, due to the slower decay time of the system with respect to the 



pump time. In all cases we again observe that the exciton number 
never exceeds 0.5 due to the same reasons as the long pump case. 

In Figure 6 we show the PL spectrum for short pump pulses. We 
again see that at the highest densities the dominant peak is around 
the cavity photon energy. As the system decays away, there is small 
shift of the peak in the direction of the LP or UP energy. At the lowest 
densities, only the LP or UP peak is visible. Apart from the decay 
dynamics, the PL signatures are fairly similar for both the long and 
short pump cases. In both cases we note that despite the high average 
densities achieved as seen in Figure 6, there is always a resonance 
peak visible at the LP and UP energies. According to the simplistic 
treatment in equation (8) we would expect that these peaks would 
only be seen at the lowest densities. We see that in fact the dynamics 
is more complicated than the simplistic picture by again considering 
Figures 2 and 3. If it were purely the case of (8) we would only see the 
density matrix elements corresponding to the highly excited states 
and zero population of the low excitation states such as Pgi,gi and 
Pxo,xo- In fact we see that for both the long and short pump-pulse 
cases, there is no reordering of the states, so that low excitation states 
are always more probable than the higher- excitation states. Thus 
there is always a strong contribution of the transitions |±, 1) 
|0) which are simply the LP and UP transitions. For this reason, 
for the typical pumping strengths that we consider, the PL spectrum 
consists of a mixture of both high and low excited states. Only the 
parts of the PL spectrum in the vicinity of the cavity photon energy 
correspond to states in the QSC regime. 

Photon PL and exciton pumping. We now discuss the case of exciton 
pumping Px with photon decay, where we fix the decay rates to k: = 1 
and y = 0. The time evolution of the PL is shown in Figure 4(b) and 
the populations in Figure 5(b). In comparison to the photon pump- 
ing case, the system requires high pump powers to reach the high 
density due to the spin operator suffering from phase-space fill- 
ing, ((7^)"|0) = 0 (n > 2). In order for the system to reach high 
densities for this pumping configuration, the pump must first create 
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an exciton in the system, then this must quickly become transferee! 
to a photon via the couphng between the cavity and the quantum dot. 
If the pump can create another exciton in the system before the 
first photon decays, then high densities may be reached. The con- 
dition required for reaching the QSC regime is therefore Px, g~>K. 
Apart from this difference, we see a quaHtatively similar spectrum to 
the photon PL/photon pumping case at high density, with the PL 
emerging most strongly at the cavity photon resonance. As the sys- 
tem decays, there is a shift of this peak towards the LP and UP 
resonances. 

For the short pump pulse case, we see an almost identical spectrum 
to the long pump-pulse case [Figures 6(b) and 7(b)]. We can there- 
fore conclude that despite the more inefficient pumping method of 
the exciton pumping, as long as the criterion Px, g~>K is satisfied, the 
PL and population characteristics are independent of the particular 
pumping scheme that is used. 

Exciton PL and photon pumping. We now turn to the excitonic PL 
and photon pumping case and fix the decay rate k: = 0 and y = 1. As 
seen in Figure 4(c), we observe a qualitatively different spectrum, 
with stronger side peaks at the LP and UP resonances, as well as the 
cavity photon resonance. We understand this by considering our 
simplified picture again. For strong photon pumping, we start once 

more in the state defined in equation (7), (a^)"|G) |0) oc 

(I + , n) — I — , n)). As the exciton PL corresponds to the loss of an 
exciton, we consider the transitions induced by operating to the 
above states 

a-\±,n)=^-{\ + ,n-\)-\-,n-\)) (10) 

Unlike the photon PL case in equation (8) where | ± , only transi- 
tions to one state, here there are two possible states to jump to. In this 
case all the transitions between \ ±,n) ^ \ — ,n — 1) have the same 
amplitude. Therefore, the most distinctive feature here is that we 
have side peaks at AE—^co + 2gy/n, in addition to A£ ~ fico for high 
density. The kind of transitions shown above are familiar from res- 
onance fluorescence, where in the high-excitation case the photolu- 
minescence spectrum reduces to the Mollow's triplet spectrum^'^^. 
The difference here is that the densities are not quite so high, such 
that we may treat the excitation field as a classical field, and also 
the non-equilibrium nature of the pumping profile. Nonetheless, 
similar spectral characteristics may be seen for the high-excitation 
regime we consider here. For continuous and high-pumping rates we 
have confirmed that our numerics reproduces the Mollow's triplet 
spectrum. 

For short pump pulses, as shown in Figures 6(c) and 7(c), the 
numerical simulations are less stable with respect to the pumping 
dynamics, hence we were only able to achieve relatively low densities 
(although still in the QSC regime) in comparison to the long pump- 
pulse cases. Nevertheless, we still see the characteristic Mollow's 
triplet spectrum evolving towards the LP and UP resonances. 
From this we conclude that both the short and long pumping profiles 
should be effective in reaching the QSC regime, giving similar spec- 
tral characteristics. 

Exciton PL and exciton pumping. We finally consider the case where 
excitons are pumped and the exciton PL is measured with y= I and k 
= 0. The time dependence of the excitonic PL is shown in 
Figures 4(d) and 6(d). We see a qualitatively similar PL spectrum 
to the excitonic PL and cavity photon pumping case, with prominent 
side peaks on either side of the cavity photon resonance. As was the 
case with photonic PL with excitonic pumping, the system requires 
high pump powers in order to reach the high-density regime, and 
requires Px, g^Y- In a similar way to the exciton PL and photon 
pumping case, all the four transitions in equation (10) occur, giving 
rise to more prominent side peaks. We thus conclude that the par- 
ticular form of the pumping from either short or long pulses, 



photonic or excitonic pumping gives qualitatively similar PL spectra 
for the excitonic PL cases. 

Discussion 

We have obtained the PL spectrum of a single quantum dot in a cavity 
pumped via injection of photons or excitons in a pulsed Gaussian 
profile. The particles decay either with the cavity photon decay rate k 
or with the exciton decay rate y depending on how we measure the 
PL, i.e. photonic PL or excitonic PL. We first find that with any choice 
of pumping scheme the quantum strong-coupling regime can be 
reached for sufficiently strong pumping strengths. While this is intui- 
tively clear for the photonic pumping case, this is also true for the 
exciton pumping, which suffers from the Pauli exclusion principle 
excluding it from the simultaneous injection of more than one exciton 
at a time. Despite this, the high density regime may be reached as long 
as the time required for conversion of excitons to photons via the 
coupling is faster than the decay rate P, g^^K, y. We also generally 
find that the PL spectrum that is measured is fairly insensitive to the 
particular pumping method used to excite the system. Our considera- 
tions of exciton and photon pumping, using both short and long 
profiles gives fairly consistent spectra throughout. Naturally, if the 
decay dynamics itself needs to be probed, then using a short pulse 
is really the only option. However, if the only goal is to reach the QSC 
regime, then either case can be used. 

The general trend of the photonic PL is that there is a peak at the 
cavity photon mode at high densities and as the system decays away it 
moves towards the UP and LP resonances in the low-density regime. 
For excitonic PL, the Mollow-triplet-like spectrum occurs as the 
density is increased. These properties agree with the theoretical pre- 
dictions of (Ref 14) at high density. Some differences to the PL 
spectrum is seen at intermediate densities (Figure 3(b) of Ref 14). 
In contrast to the many-exciton model, where the photolumines- 
cence smoothly transitions from the LP to cavity photon energy, in 
the single quantum dot case the PL discontinuously shifts between 
the two limits. We have confirmed that by increasing the number of 
excitons in the system the transition gradually becomes smoother. 
The discontinuous spectrum can therefore be seen as a signature of 
the quantum dot system, although in a typical experimental situation 
it is more likely that many quantum dots are excited simultaneously 
in order to have sufficient signal amplitude in the measurement. 

Methods 

Time evolution of the density matrix. Due to the open-dissipative nature of the 
quantum dot-microcavity system, the PL characteristics are modeled by a master 
equation following Refs. 9-11. We consider the reduced density matrix 

j^p=^[p, TCsys] + I {2apa' -a'ap- pa' a) 

+ -{2(T~ p(T^ —o^o~p — po^o~) 

H — {2a'pa — aa'p — paa') 

+ ?^(2G+pa- -a-G+p-pG-G+), 
where we assume that the pump has a Gaussian time profile 

Pc(0,f'x(0 = f'o exp(-^^^). (12) 

Using the basis | in) where / = G is the state with no exciton and / = X is the state with 
an exciton, we may define the matrix elements of the reduced density matrix 

Pin,im = {Mp\j^)^ (13) 

where n, m is the photon number. The time evolution of the density matrix elements 
are 

JfPGn,Gn = ^gVn{pGn,Xn-l - PXn-l,Gn) +yPxn,Xn 

- K [npGnfin " + 1)Pg«+ l,Gn+ l] ( H) 

-PxpGn,Gn+PcnpGn- l,Gw-l 

-Pc(«+l)PGn,G«' 
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I i , — , , 

-PXn,Xn = ^gVn^l(Pxn,Gn+l-PGn+l,Xn) -yPXn,Xn 

-^[npXn,Xn - («+ 

+ PxPGn,Gn+PcnPxn- 

-Pc{n+l)Pxn,Xn^ 



J^PGn,Xn-\ = ^ [gVn{pGn,Gn- Pxn-l,Xn-l) +^PGn,)i 



(15) 



-[{y + K{2n-l)+Px)l2]pcn,, 



+ ^^/n(n + l)pGn+\,x 



(16) 



dt 



+ PcVn(n^PGn-i,Xn-2 
-y (2«+l)pG„,Xn-l. 

Pxn-l,Gn=- ^ [gVn{pGn,Gn- Pxn-l,Xn-l) 
+ ^PXn-l,Gn] 

-[iy^Ki2n-l)^Px)/2]pxn-i,Gn 



(17) 



'Xn-2,Gn-l 



+ K^/n{n+l)pxn,Gn + 

+ Pc\/n{n-l)p_ 
-^(2n+l)pxn-i,Gn- 

Imposing a maximum photon number n^^^ we may evolve this set of equations 
which form a set of closed ~ 4/imax equations. In order to investigate the system 
dynamics, we truncate up to finite times instead of calculating up to the stationary 
limit. 

Two-time correlation function. In the quantum dot system, two types of 
photoluminescence mechanism can be considered: the particles are supplied from the 
pump laser and they decay due to leakage through the microcavity mirrors, or via 
exciton recombination^^. We call these processes the photonic and excitonic PL, 
respectively. The PL is given by^ 



S(r,v) = 



-US 

Jo 



dx e'"'G^'\t,T), 



(18) 



where G^^^f, t) is the two-time correlation function, which for the quantum dot 
system, is either 



G(.'>(t,T) = (flt(t + t)fl(t)), 
G<;"(t,T) = (a+(t + T)a-(0), 



(19) 
(20) 



depending on which type of PL is being calculated. The Fourier transformation of 
G^^\t,T) gives the photonic PL which is emitted by the leakage from the imperfect 
mirror, while thatof G^^V.-^) 

gives the excitonic PL which is emitted due to the 

recombination of excitons. 

Following Refs. 9-11, in the interaction picture we may conveniently define the 
following operators 

4„W = |Gn + l>(G«|e'('"»-^W«, 
4„(T) = |X« + l>(Xn|e'(''»-*W», 
(T+{r) = \Xn){Gn\e"'"-"''>, 
C„(T) = |G>i + l>(Xn-l|c'<'"»-''^>'/*. 



(21) 



In terms of the new operators defined in equation (21), we can rewrite the two-time 
correlation function as 

G^^\t,T)=Y,V^{{4„it + rHt)) + {al{t + x)a{t))), (22) 



(23) 



In order to calculate the two-time correlation function, we need to obtain the 
expectation values of two operators at different times. The time evolution of the 
expectation values for each operators is then 



^ (4(^)> = (4(^)> (^«- A) - ^ (2. + 1) -Px 
+ (4„+iW)'<^\/(" + l)(" + 2) 

-^[{4n-2i^))fPcVn{n-l) 

+ (^:+i(t))k:(«+ 1) - {al^^ M)^gV^ 
+ {<-iM)Pcn-{(J^ (t))Pc(«+ 1), 



(24) 



(25) 



4n-A^)) = {4n-,M)Px-{<iT))^gV^ 

+ (4„-iW) ^(^«-A)-y-^(2«-l) 

+ {4nMyVn{n + l) + (?„(t)> ^gV^ 
+ {al.2M)PcVn{n-l) 

-(4„-iW>Y(2« + 1), 

^ (QnM) = - (4« W) ^^^+ {4n-li^))^gV^ 

+ (g„(T)>t(/za3-A)-^^^^-K:^ 



(26) 



+ (?„+i(t))/cv'«(« + 2) 



(27) 



+ {Qn-iM)PcVin-l){n+l) 
-{Q„{T))Pc{n + l). 

Applying the quantum regression formula to equations (24)-(27), we can compute 
and obtain the following sets of two-time correlation functions^^. For the calculation 
of G^ (t,T), we need to obtain 



{{al{t + T)a{t)),{a:{t + T)a{t)), 

(4„-i(^ + ^)«(0),fen(f + T)a(0>} 

and, for the calculation of G^^^ (t,T), 

{[al{t + r)a-{t)),{a+{t + r)a-{t)), 

{al_,{t + T)^-{t)),{Ut + T)^-{t))] 

is needed. The above equations satisfy the initial conditions 

(4„W«W) = V^Pg„+i,g„+i(0, 
((rt(<)«(0) = \/^PG„+i,x„We'^'''*, 
{<i„-i(tHt)) = V^Px„.x„(t), 

(?„(t)«(<)) = \/Spx„,G„+i(f)e-"'/", 



(28) 



(29) 



(30) 



and 



(aL(0^(0>=Px.,G„+l(Oe-'^^/^ 

{<jlitHt))=Pxn.Xnit), 

(aL_i(O^W)=0, 
{g„{tMt))=0, 



(31) 
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for the calculations of g[^^^ {t, t) and G^^^ {t, t), respectively. Here, t is calculated up to 
the stationary limit so that we can compute the two-time correlation functions. 
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